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METHOD AND SYSTEM FOR ESTIMATING INTERCONNECT DELAY 



TECHNICAL FIELD OF THE INVENTION 

This invention relates generally to the field of 
system design and more specifically to a method and 
system for estimating interconnect delay. 
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BACKGROUND OF THE INVENTION 

Accurate calculation of interconnect delay is 
critical to the design of high speed integrated circuits. 
Signals travel along paths that include the interconnects 
5 and gates of an integrated circuit . The interconnects 

contribute to the delay of the signal along the path. 
Accurate calculation of interconnect delays, however, is 
typically time-consuming and expensive. Accordingly, 
calculating interconnect delay has posed a challenge for 
10 integrated circuit designers. 
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SUMMARY OF THE INVENTION 

In accordance with the present invention, a method 
and system for estimating interconnect delay are provided 
that substantially eliminate or reduce the disadvantages 
5 and problems associated with previously developed systems 

and methods . 

According to one embodiment of the present 
invention, a method for estimating interconnect delay is 
disclosed that includes determining inductance of an 

10 interconnect. A transfer function is determined using the 

inductance, and two poles of the transfer function are 
determined. An interconnect response is estimated using 
the two poles, and an interconnect delay is estimated 
from the interconnect response. 

15 According to one embodiment of the present 

invention, a system for estimating interconnect delay is 
disclosed that includes a memory that stores information 
about an interconnect . A processor determines an 

inductance of the interconnect from the information, a 

20 transfer function using the inductance, and two poles of 

the transfer function. The processor also estimates an 
interconnect response using the two poles, and an 
interconnect delay from the interconnect response 

Embodiments of the invention may provide numerous 

25 technical advantages. A technical advantage of one 

embodiment is that the interconnect delay is determined 
using the inductance of the interconnect . Global 
interconnects typically have large cross- sect ions and are 
driven by drivers with small resistance. This may result 

30 in an inductive impedance that is greater than a 

resistive impedance and that affects interconnect delay. 
Accordingly, interconnect delay is estimated using 
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inductance in order to improve the accuracy of the 
estimation. 

Another technical advantage of one embodiment is 
that the interconnect delay is estimated according to a 
damping response of the interconnect. When inductance is 
taken into account, the damping response may also be 
determined, which provides a better description of the 
interconnect. Accordingly, estimating interconnect delay 
according to the damping response provides an improved 
estimate . 

Another technical advantage of one embodiment is 
that the embodiment may be used to estimate interconnect 
delay for system design, interconnect design, circuit 
design, board design, and backplane design, among other 
areas of computer design. 

Other technical advantages are readily apparent to 
one skilled in the art from the following figures, 
descriptions, and claims. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

For a more complete understanding of the present 
invention and for further features and advantages, 
reference is now made to the following description, taken 
5 in conjunction with the accompanying drawings, in which: 

FIGURE 1 illustrates a routing tree with 
interconnects for which interconnect delay may be 
estimated; 

FIGURE 2 is a block diagram of a system for 
10 estimating interconnect delay; and 

FIGURE 3 is a flowchart of a method for estimating 
interconnect delay . 
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DETAILED DESCRIPTION OF THE DRAWINGS 

Embodiments of the present invention and its 
advantages are best understood by referring to FIGURES 1 
through 3 of the drawings, like numerals being used for 
5 like and corresponding parts of the various drawings. 

FIGURE 1 illustrates a routing tree 10 with 
interconnects 102 for which interconnect delay may be 
estimated. Routing tree 10 represents, for example, the 
layout of an integrated circuit such as a large-scale 

10 integration (LSI) , very-large-scale integration (VLSI) , 

super- large -scale integration (SLSI) , or ultra-large- 
scale integration (ULSI) integrated circuit. Routing 
tree 100 includes nodes 106 and interconnects 102, 104, 
and 109 coupling nodes 106. Nodes 106 represent elements 

15 such as transistors of an integrated circuit, and 

interconnects 102, 104, and 109 represent connections 
such as wires between the elements. Interconnect 104 
from node S 106a to node L 106h, where node S 106a is a 
source node and node L 106h is a load node, includes 

20 interconnects 102 and nodes 106. 

Response describes the behavior of an interconnect 
in accordance with an applied signal, for example, the 
response of interconnect 102d describes the behavior at 
node B 106e in response to a signal applied at node A 

25 106d. Delay describes the time interval between the 

instant a signal is applied to an interconnect and the 
instant the interconnect acts accordingly, for example, 
the time interval between the application of a signal at 
node A 106d and a response to the signal at node B 106e. 

30 A transfer function of an interconnect is used to derive 

the response and delay of the interconnect . Parameters 
of an interconnect may be used to determine attributes of 
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the interconnect, which in turn may be used to derive 
coefficients, or alternatively moments, of a transfer 
function. Parameters of an interconnect such as a wire, 
include the length h, width w, thickness j, and 
5 dielectric coefficient of the wire. The distance d 

between interconnects may also be a parameter. 
Attributes may include the resistance R, capacitance C, 
and inductance L.of the interconnect. 

FIGURE 2 is a block diagram of a system 100 for 

10 estimating interconnect delay. System 100 receives 

parameters of an interconnect 102, for example, a length 
h, width w f and thickness j of interconnect 102 and a 
distance d between interconnects, and estimates an 
interconnect delay T for interconnect 102. 

15 System 100 includes a processor 110 coupled to a 

memory 111, an inductance module 112, a transfer function 
module 114, an overdamped module 120a, an underdamped 
module 120b, and a critically damped module 120c. The 
operations of the processor and modules of system 100 may 

20 be performed by any number, combination, or arrangement 

of processors or modules, and may be performed by any 
combination of hardware, software, or both. Processor 110 
manages system 100 to estimate interconnect delay and 
stores data in memory 111. Inductance module 112 

25 receives parameter information about interconnect 102 

from processor 110, and computes an inductance L of 
interconnect 102. Transfer function module 114 

determines a transfer function H(s) using inductance L 
calculated by inductance module 112, and other 

30 interconnect attributes such as resistance * R and 

capacitance C calculated by processor 110. Processor 110 
determines the poles s 2 and s 2 of transfer function H(s) , 
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and determines a damping case, for example, an overdamped 
response, an underdamped response, or a critically damped 
response, from poles s 2 and s 2 . 

Overdamped module 12 0a receives poles s 2 and s 2 if 
5 the damping case is an overdamped response. A response 

module 122a computes an overdamped response v(t), and a 
delay module 124a computes a interconnect delay T. 
Underdamped module 12 0b receives poles s x and s 2 if the 
damping case is an underdamped response. A response 
10 module 122b computes an underdamped response v(t), and a 

delay module 124b computes a interconnect delay T. 
Critically damped module 120b receives poles s 2 and s 2 if 
the damping case is a critically damped response. A 
response module 122c computes a critically damped 
15 response vftj, and a delay module 124c computes a 

interconnect delay T. 

Unlike known systems, system 100 for calculating 
interconnect delay takes interconnect inductance into 
account. Additionally, system 100 calculates 

20 interconnect delay based on a damping response, unlike 

known systems. Consequently, system 100 provides 

accurate estimates of interconnect delay . 

FIGURE 3 is a flowchart of a method for estimating 
interconnect delay. The method begins at step 210, where 
25 interconnect parameters are received by processor 110. 

Interconnect parameters include length h, width w, and 
thickness j of interconnects 102, and a distance d 
between interconnects 102 . 

Interconnect inductance, such as the self and mutual 
30 inductance of the interconnects, is calculated by 

inductance module 112 at step 212. In one embodiment, 
interconnects 102 may include signal and shield lines. 
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The width w of the signal line is approximately equal to 
the spacing between adjacent signal lines, and the width 
of a shield line is approximately twice the width w of a 
signal line. Self -inductance L se if for the signal and 
5 shield .^lines may be computed from interconnect parameters 

length h, width w, and thickness j using Equation (1) : 

L eIf =-^-*h*[0.5 + ln(— + 0.2235(w + j) ] (1) 

2 * 71 w + J 

The mutual inductance M 12 between two lines separated by 
a distance d is given by Equation (2) : 

10 M = ^° * h * [ln(— ) + - 1 + — ] (2) 

2 * 7i d h 

The inductance L Ioop of a current loop between a signal 

line and a shield line may be computed using Equation 
(3) : 

■^loop = ^self + ^self ~~ ^ sig2shld ~ ^ shld2shld ( 3 ) 

15 L"^ represents the self -inductance of a signal line, and 

L s s h e \ d f represents the self -inductance of a shield line. 

M S ig 2 shid represents the mutual inductance between the 
signal and shield lines on both sides of the signal line, 
and M S hid23hid represents the mutual inductance between two 
20 shield lines on both sides of the signal line. The 

inductance L of an interconnect is the sum of the self- 
inductance and the loop inductance. 

A transfer function H(s) is calculated using 
inductance L at step 214. Transfer function module 114 
25 provides processor 110 with a transfer function H(s) , 

which is expressed by Equation (4) : 

1 



H (s) = 



1 + b x s + jb 2 s 2 (4 ) 

= 1 + sM, + SM 
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The coefficients jb 2 and b 2 are functions of interconnect 
102 and receiving node 106 attributes, as expressed by- 
Equations (5) : 

RC 

1 = T" + * c " 

{RC) 2 r 2 cc„ LC (5) 

b = + - + + LC N 

24 6 2 

5 where C N is the receiving node 106 capacitance, and R is 

resistance, C is capacitance, and L is the inductance of 
the interconnect 102. Inductance L is received from 
inductance module 112, and resistance R and capacitance C 
and C N are computed by processor 110 from received 
10 interconnect parameters. 

Poles Sjl and s 2 of the transfer function H(s) are 
determined at step 216. Processor 110 receives transfer 
function H(s) and determines poles s 2 and s 2 of the 
transfer function H(s), as expressed by Equation (6): 



2 - b ± Jb 2 - 4b 

15 s l2 = , 2 = 2Li L (6 ) 

- AT X ± A /4M 2 - 3M' 2b 2 



A damping case is determined from poles s 2 and s 2 at 
step 218. If the poles are real, the damping case is 
determined to be an overdamped response. The condition 
for the poles to be real is given by Equation (7) : 

20 A /4M 2 - 3Ml = - 4Jb 2 > 0 (7) 

If the poles are complex, the damping case is determined 
to be an underdamped response. The condition for the 
poles to be complex is given by Equation (8) : 

^bl - 4b 2 < 0 (8) 

25 If there are double poles, the damping case is determined 

to be a critically damped response. The condition for 
double poles is given by Equation (9) : 
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^Jb 2 - 4jb 2 = 0 (9) 

Overdamped module 120a receives poles s 2 and s 2 at 
step 22 0, if the damping case is determined to be an 
overdamped response. A response is estimated at step 
5 222 . Response module 122a receives the poles and 

computes an overdamped response v(t) using Equation (10): 

v(t) = V 0 (1 ^ e SlC + — ^ — e fllC ) (10) 

where V 0 is the steady state voltage. A delay is 
estimated from the overdamped response v(t) at step 224 
10 by delay module 124a. The contribution of s 2 to 

overdamped response v(t) may be neglected, since s 2 is 
smaller than Si . Delay module 124a uses Equation (11) to 
compute the delay T: 

T = In 1 (1 — v t J * J- = in 1 (1 - v th ) * f** (11) 

Nil b, - ^bl -4Jb 2 

15 where v th is the threshold voltage for the damping 

response . 

An underdamped module 12 0b receives poles s ± and s 2 
at step 230, if the damping case is determined to be an 
underdamped response. A response is estimated at step 
20 232. Response module 122b computes an underdamped 

response using Equation (12): 

f / _.2 . o 2 \ 



v(t) = v 



Jet 2 + (3 

1 - * — — — e" at * sin ((3 1 + p) 



(12) 



M n J3M - 4M, B 

where a = — * , p = ^ 1 2 - , and p - tan" 1 (-^) . 

2(M 2 - M 2 ) 2(M 2 - M a ) a 

The value of a sine function sin(pt + p) is calculated at 
25 step 234 to determine whether a sine delay function or an 

exponential delay function is to be use to estimate the 
interconnect delay . 
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If sin(pt + p) is zero or greater at step 234, the 
method proceeds to step 2 36, where a sine function is 
used to estimate the delay. Delay module 124b computes a 
delay T from response v(t) by solving response v(t) 
5 recursively for the time variable t. A next delay T nex t is 

calculated using a current delay value T cur for the time 
variable t of a sine function, as expressed by Equation 
(13) : 



10 



15 



20 



T = — * In 
a 



* sin (J3 * T cur + p ) * 



(13) 



a - v t j 

In the initial delay T 2 calculation, an Elmore delay value 
T ED may be used as the current delay value for the time 
variable t of the sine function, as expressed by Equation 
(14) : 



T x = — * In 
a 



_(1 ~ v 

where the Elmore delay 



(14) 



* sin (p * T ED + p) * V 1 + a 7P 2 

* jb x , which is 



Ted = In 



(1 " V t J 

7^=0.22^ for a 20% threshold as the starting time. The 
method then proceeds to step 23 8 to determine whether the 
estimate is acceptable. 

Referring back to step 234, if the sine function 
sin(Pt+p) is negative, the method proceeds to step 237, 
where the delay is estimated using an exponential 
function, as expressed by Equation (15) : 



T — — — =z K * 



2£>, 



(15) 
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d - 

where K = = p. A negative sine function 

F¥ 

sin(pt+p) indicates that the interconnect inductance has 
a stronger effect than the interconnect resistance. The 
method then proceeds to step 23 8 to determine whether the 
5 estimate is acceptable. 

At step 238, if the difference is not within the 
predetermined range, for example, a range where T next e 
[0.5T cur/ T cur ] , then next delay T nex t is not an acceptable 
estimate for the interconnect delay, and the method 
10 returns to step 234 to determine whether to use the sine 

delay function or the exponential delay function. 

If sin(Pt + p) is zero or greater at step 234, the 
method proceeds to step 236, where the sine function is 
used to estimate the delay. At step 236, another next 
15 del ay T next is computed using the previously computed delay 

for the value of the time variable of the sine function, 
as expressed by Equation (16) : 



1 

T 2 = — * In 
a 



* sin (p * T x + p) * Jl + a 2 /p 2 

(1 " v th ) 



(16) 



The method then proceeds to step 23 8 to determine * whether 
20 the estimate is acceptable. 

If the sine function sin(PTx+p) is negative at step 

234, the method proceeds to step 237, where an 

exponential function is used to estimate the delay. 

Equation (15) , where the previously computed delay is 
25 used for the value of the time variable of the 

exponential function, that is, where K = Vth ^ e p t is 
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used to compute the delay. The method then proceeds to 
step 238 to determine whether the estimate is acceptable . 

At 238, if the difference between the next delay T next 
and the current delay T cur is within a predetermined 
range, then next delay T next is an acceptable estimate for 
the interconnect delay at step 238, and the method 
terminates . 

A critically damped module 120c receives poles s 2 and 
s 2 at step 240, if the damping case is determined to be a 
critically damped response. A response v(t) is estimated 
at step 242 . Response module 122c computes a critically 
damped response using Equation (17) : 



v(t) = V 



f 2t A 

1 - e ts ' - — e ts * 

V h r J 



(17) 



A delay T is estimated at step 244 by delay module 
by 124c. Delay T is calculated from response v(t) by 
solving response v(t) recursively for the time variable 
t. A next delay T next is computed using a current delay 
value T cur for the time variable t of an exponential 
function, as expressed by Equation (18) : 

2b ? __ . Jb, 



T — K * ~ K * — (IS) 

^next A next , ^ next _ V -L o / 



'1 



where K next = In 



2T 
(1 + — 



1 - V tn ^1 



In the initial delay T 2 calculation, an Elmore delay 
T ED is used as the current delay value for the time 
variable t of an exponential function, as expressed by 
Equation (19) : 

2b, jfcx 

T i = K i * — ~ = K i * — (19) 
*>i 2 



S^^ZKET J^AI 

n /1 mi r\r\\ 
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where K = In 



2T 
(1 + — =-) 



For a 20% threshold, the 



.(1 " vj 

Elmore delay is T ED = 0.2 2*bi. 

If the difference is not within a predetermined 
range, for example, a range where T next e [0.5T cur , T cur ] , 
the method returns to step 244, where another delay is 
computed. The next delay is determined using the 

previously computed delay as the time value for the 



exponential function, that is, K 2 = In 



* (1 + — L ) 



.(1 - v^) b, 

If the difference between the next delay T next and the 

10 current delay T cur is within a predetermined range, next 

del Tnext is an acceptable estimate for interconnect 
delay at step 238, and the method terminates. 

A technical advantage of one embodiment is that the 
interconnect delay is determined using the inductance of 

15 the interconnect. Global interconnects typically have 

large cross- sect ions and are driven by drivers with small 
resistance. This may result in an inductive impedance 
that is greater than a resistive impedance and that 
affects interconnect delay. Accordingly, interconnect 

20 delay is estimated using inductance in order to improve 

the accuracy of the estimation. 

Another technical advantage of one embodiment is 
that the interconnect delay is estimated according to a 
damping response of the interconnect. When inductive 

25 impedance is taken into account, the damping response may 

also be determined, which provides a better description 
of the interconnect. Accordingly, estimating interconnect 
delay according to the damping response provides an 
improved estimation. 
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Although an embodiment of the invention and its 
advantages are described in detail, a person skilled in 
the art could make various alterations, additions, and 
omissions without departing from the spirit and scope of 
the present invention as defined by the appended claims. 



